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Abstract 

We study the exclusive decay of 1 heavy quarkonium into one photon and two pions 
in the kinematic region, where the two-pion system has a invariant mass which is much 
smaher than the mass of heavy quarkonium. Neglecting effects suppressed by the inverse of 
the heavy quark mass, the decay amplitude can be factorized, in which the nonperturbative 
effect related to heavy quarkonium is represented by a non-relativistic QCD matrix element, 
and that related to the two pions is represented by a distribution amplitude of two gluons 
in the isoscalar pion pair. By taking the asymptotic form for the distribution amplitude 
and by using chiral perturbative theory we are able to obtain numerical predictions for the 
decay. Numerical results show that the decay of J/ij: can be observed at BEPC and at 
CESR. Experiment observation of this process in this kinematic region at BEPC and CESR 
can provide information about how gluons are converted into the two pions and may supply 
a unique approach to study / = s-wave vrvr scattering. 
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Recently it has been proposed to study production of two pions in exclusive processes like 
7 + 7* ^ TT + TT 1^, ^ and like 7* + /;,^/;, + 7r + 7r0, ^. These processes will enable us to study 
how quarks and gluons, which are fundamental dynamical freedoms of QCD, are transmitted 
into the two pions. In the kinematic region, where the two pions have a much smaller invariant 
mass than the virtuality of the virtual photon, the scattering amplitudes of the processes take a 
factorized form. In this factorized form the transition of partons into the two pions are described 
by distribution amplitudes of partons in the two pions, which are defined with twist-2 operators. 
For 7 + 7* —>■ TT + 77, at the tree-level, only the distribution amplitude of quark appears in 
the scattering amplitude, the distribution amplitude of gluon appears at loop levels or through 
evolution of distribution amplitudes[^ |^, |^, while for 7* + /i — > /i + vr + vr, at the tree-level, the 
distribution amplitude of quark as well as that of gluon contribute to the scattering amplitude, 
and the produced charged pion pair is in isospin / = or / = 1 states]^, In all of the above 
processes, the produced system of two charged pions will be dominantly in an isospin / 7^ state. 
This may make it difficult to extract the distribution amplitude of gluon from experimental data, 
because the two pions produced through gluon conversion are in a / = state. In this letter we 
propose to study the radiative decay of J/ip into two pions, where the distribution amplitude of 
gluon appears at the tree-level and the produced two pions are dominantly in a / = state. This 
makes the extraction of gluon content of a two-pion system relatively easier in experiment. Beside 
this the decay also provides a possibility to study / = s-wave vrvr scattering. The decay can be 
studied with the data sample of 5 x 10^ J/ip^s whose collection will be completed with BES at the 
end of this year. 

We consider the decay in the kinematic region where the two pions have a invariant mass 
which is much smaller than the mass of J/ip. In the limit of rric —>■ 00 the decay amplitude can 
be factorized, where the nonperturbative effect of J/ip can be presented by a NRQCD matrix 
element, while the nonperturbative effect of the two pions is represented by the same distribution 
amplitude of gluon appearing in the exclusive processes mentioned above. Corrections to this limit 
may systematically be added. With the model for the distribution amplitude of gluon, developed 
in 0, 0, we obtain numerical predictions for the decay in the heavy mass limit, and they indicate 
that the decay is observable at BES. 

We study the exclusive decay in the rest frame of J/ip: 



The momenta are indicated in the brackets. We denote k = k+ + k_ and m^^ = k'^. At leading 
order of QED, the S-matrix element for the decay is 



where Qc is the electric charge of c-quark in unit of e, c{z) is the Dirac field for c-quark. At 
leading order of QCD, two gluons are emitted by the c- or c-quark, and these two gluons will be 
transmitted into the two pions. Using Wick theorem we obtain: 



J/^(P) ^ 7(A;^) + 7r+(/c+) + 7i-{k_). 



(1) 



{f\S\t) = -iQ, e 1 <fzlaT^^Ti~\A^{z)-c{z)^>'c{z)\Jl^) 



(2) 



% 1 




,i{k.,-z+k2-y) 



2 2 



1 



X (0|9(a;i)Q(yi)|JM (7r+7r-|G;ja;)Gt^(0)|0) 

X [5\x - x^)5\z - y{)i^^SF{x - yh'^Spiy - z)Y + ■■■], (3) 

where /c2 is the momentum of one of emitted gluons, Sf{x — y) is the Feynman propagator of 
c-quark, the dots in the square bracket denotes other five terms. In the hmit of rric oo, a c- 
or c-quark moves with a small velocity v, this fact enables us to describe nonperturbative effect 
related to J/ip by NRQCD, in which a systematic expansion in v is employed. This expansion 
is extensively used for inclusive decays 0. This expansion can also be used here for the matrix 
{0\cj{xi)ci{yi)\J/ip) , in which we expand the Dirac field c{x) with corresponding NRQCD fields. 
With the expansion we obtain: 

{0\cJ{x,)c.{y,)\J/^P) = - J (P+ V P_)., (0|xVV| W) e-^^-(^^+^^) + 0{v'), (4) 

o 

where xKi^) is the NRQCD field for c(c) quark, a\l = 1,2,3) is the Pauli matrix, and 

p = K, 0,0,0). (5) 

The matrix {0\x^a'"ip\J/ip) is proportional to the polarization vector e\J/ip) at the considered 
order. In this paper, we neglect the contribution from higher orders in w, the momentum of J /if) 
is then approximated by 2p. It should be noted that the effect of higher order in v may be added. 
Using Eq. @) we can write the S-matrix element as 

{f\S\i) = ^Qc e gl {27cr 5\2p -k-k,) e;(7) (0|xVV| W) 

/ ^,H'''''''{p,k,k,)T,,,,{k,k,), (6) 

and 

T'"'{k,ki) = J d^x e"^^^^ (7r+7r-|G"^'^(a;)G''''^(0)|0), (7) 

in which if'^'^'i^i^p^ k^ ki) is a perturbative function, which is the amplitude for a cc pair emitting a 
photon and two gluons. The contributions in Eq. can be represented by Feyman diagrams. One 
of them is given in Fig.|I|, where the kinematic variables are also indicated. The nonperturbative 
object T^^^{k, ki) describes how two gluons are converted into 2tt. 

In the heavy quark limit the two-pion system will have a large momentum fc, we take a 
coordinate system in which the direction of k is the 2;-direction. For convenience we will work in 
the light-cone coordinate system, in which the components of k are given by 

k^" = {k+,k~,Q), k+ = {k^ + k^)/V2, k- = {k'^ -e)/V2. (8) 
In the light-cone coordinate system we introduce two light cone vectors and a tensor: 

= (0,1,0,0), P = (1,0,0,0), 
4^ = gf^^ -n^'r -n'^l'', (9) 
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Figure 1: One of the Feynman diagrams for the exclusive decay of J/ip into one photon and two 
pions. 



and we take the gauge 



n 



G{x) 



0. 



(10) 



The x-dependence of the matrix element in Eq.(7) is controlled by different scales. The x~- 
dependence is controlled by k~^, which is large, while the x'^— and XT-dependence are controlled 
by the scale k~ and Aqcd, which are small in comparison with k'^. With this observation we 
can expand the matrix element in x~ and in xt- At the leading order only twist-2 operators 
contributes to the matrix element. We will neglect higher orders in the expansion, i.e., we only 
keep contributions of twist-2 operators. Then we obtain: 



r^^''ik,k,) = i27ry SikT) S\k,T) 



1 



1 



k+ xi{l — xi) 



1 



(11) 



with 



$ (a;i,C,m^^) 

Xi 



1 

k+ 
ki ■ n 

k ■ n 



dx 

c 



X (7r+7r-|G'"'+^(x-n)G'%(0)|0), 



/c+ ■ n 
k ■ n 



(12) 



$*^(a;i, C, ""^TTTr) is the gluonic distribution amplitude which describe how a pion pair is produced 
by two coUinear gluons, it represents a nonperturbative effect and can only be calculated with 
nonperturbative methods or extracted from experiment. As it stands, it is gauge invariant in the 
gauge n ■ G{x) — 0. In other gauges we need to supply a Wilson line operator to make it gauge 
invariance. With Eq.(ll) it is straightforward to obtain the -S-matrix element at the tree- level in 
our approximation: 



{m^) = 



—I 
I2ml 



e gl (27r)^ 5\2p - k - fc,)e*^(7)(0|xVV| 



f 
Jo 



dxi 



1 8ml{l - xi) - ml^jl ^ 2xi] 

xi{l — xi) Aml{l — xi) — ml^{l — 2xi] 



^'^{xi,C,rn^n)- 



(13) 



In the calculation we have kept the mass 111^^^^. Usually the effect of m^7r should be combined with 
effects of twist-4 operators as a correction to the above result. We have checked that striping the 
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TTT-TrTr teims changes the results by less than 5%. In Eq.(13) the nonperturbative effect related to 
J /if) and that to the two-pion system are separated, the former is represented by a NRQCD matrix 
element, while the later is represented by the matrix element of gluonic operators, which is defined 
in Eq.(12). The NRQCD matrix element is related to the wave-function of J /ip in potential models 
and can be estimated with these models. It can also be calculated with lattice QCD or extracted 
from experiment. Little is known numerically about the matrix element of gluonic operators, 
i.e., the function C, ""^TTTr)- A detailed study of properties of C, ""^TrTr) can be found 

in [|], ^ For our numerical prediction we will use their results for $'^(a;i, C, 'tT'tttt), in which 
asymptotic form of C, J^tttt) is taken as an Ansatz for $'^(a;i, C, "^tttt)- It should be noted 

that the renormalization scale should be taken as 2mc in our case. Because it is not very large, 
the actual shape of C, ""^tttt) may look dramatically different than that of the asymptotic 

form. Keeping this in mind we take the form $'^(a;i, C, ""^tttt) as that given in 0: 

C, m^n) = -60 xl (1 - 

(14) 
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Figure 2: The differential decay branching ratio d 'BYi{J / ip) / dm-j^T^ as a function of m,r7r in unit of 
lO-^GeV"^ 

where 9^ is the scattering angle of t^^ and (3 is the velocity of 7r+ in the center of mass system of 
the two pions, they are related to C and m^jr by 



/3 cos^c = 2 C - 1, 13 



\ 



4 ml 



C is a constant and takes C = 1 + h ml + 0{m^) with b ~ — 1.7GeV ^ 
fraction carried by gluons in the pion, its asymptotic value is 



4 Cf 



4 



(15) 

, 1^, M2 is the momentum 

(16) 




M ( GeV ) 



Figure 3: The decay branching ratio of J/?/' — 7 + tt+tt as a function of M, the upper cut of 



foifnmr) and /2("^7r7r) are the Omnes functions for I = S- and D-wave tttt scattering, respectively. 
The Omnes function f2{jn-K-K) is dominated by the /2(1270) resonance resulting a peak at m^^r = 
1.275GeV, while the Omnes function /o("^7r7^) in the relevant m^^r region we studied (m^7r < 
0.70GeV) can be calculated by the chiral perturbative theory, the result is 



1 + 



+ 



2m! 



ml 



where 



192 7r2 /2 32 
93MeVis the pion decay constant. 



(3 In 



1 + 



+ 2 + mf3 



(17) 
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Figure 4: The differential decay branching ratio d BR(T)/(im^jr as a function of mj^^ in unit of 
lO-^GeV^^ 

With these results we are able to predict the shape of the differential decay branching ratio 
d BR{J / ip) / dmT,T, numerically where we use the leptonic decay of J/ip to determine the NRQCD 
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matrix element. The result is presented in Fig.2, where we take the parameters: nic = l.SGeV 
and as{2mc) = 0.31. We also calculate the decay branching ratio, where we use a upper cut M 
for mTTTT, hence the decay branching width is a function of M: 



2m^ dfri-n 



The numerical result for this decay ratio is given in Fig. 3. For M > 0.38GeV, we can read 
from Fig. 3 that the branching ratio is larger than 10~^, for M = 0.70GeV the branching ratio is 
9.3 X 10""^. Therefore it is likely that the decay mode can be measured using the data collected 
at BES. In the kinematic region we consider the dominant contribution to the differential decay 
width is from /o(^«7r7r)- The contribution from /2(?7i7r7r) is less than 10%. This gives a possibility 
to study / = s-wave tttt scattering in the decay. Similar results for T-decay are also obtained, 
which are presented in Fig.4 and in Fig. 5. For m,r7r < 0.70GeV, the branching ratio is 2.0 x 10^^. 
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Figure 5: The decay branching ratio of T — >• 7 + tt+tt^ as a function of M, the upper cut of m-^n 
in unit of 10~^. 

It should be noted that the two-pion system produced in the decay is in a S-wave state in 
our approximation. At twist-2 level, the produced system can also be in a D-wave state, the 
corresponding tensor distribution amplitude can be found in [0]. However the decay amplitude 
for a D-wave state is kinematically suppressed, because it is proportional to {k^ — kS)^{kj^ — k_Y 
and the relative momentum (/c_|_ — kJ) of the two pions is a small quantity in the kinematic region 
considered here. If one replaces the two-pions system by the tensor meson /2(1270), one also finds 
that the decay amplitude for /2 with the helicity ±2 is kinematically suppressed^]. Although the 
effect of the D-wave state is suppressed, it can lead to some nonzero azimuthal asymmetries which 
may be observable in experiment. 

To summarize: In this letter we have studied the exclusive decay of 1 heavy quarkonium 
into one photon and two pions. At leading twist, the S-matrix of the decay takes a factorized 
form, the nonperturbative effect related to heavy quarkonium and that related to the light hadrons 
can be separated, the former is represented by a non-relativistic QCD matrix element, while the 
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later is represented by the two-gluon to isoscalar two-pion distribution amplitude. By taking the 
asymptotic form for the distribution amplitude and by using chiral perturbative theory we are 
able to obtain numerical predictions for the decay. For J /ip , the branch ratio for m^rTr < O.TOGeV 
is 9.3 X 10""^, while for T, the branch ratio for m^r^ < 0.70GeV is 2.0 x 10~^. The decay of J/ip can 
be studied at BES and the study can provide a possibility to study 7 = s-wave tttt scattering. 
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